Functions of Noncommuting Operators 
in an Asymptotic Problem 
for a 2D Wave Equation with Variable Velocity 
and Localized Right-Hand Side 

Sergei Dobrokhotov, Dmitrii Minenkov, Vladimir Nazaikinskii, 
and Brunello Tirozzi 

Dedicated to Vladimir Rabinovich 



Abstract. In the present paper, we use the theory of functions of noncommut- 
ing operators, also known as noncommutative analysis (which can be viewed as 
a far-reaching generalization of pseudodifferential operator calculus), to solve 
an asymptotic problem for a partial differential equation and show how, start- 
ing from general constructions and operator formulas that seem to be rather 
abstract from the viewpoint of differential equations, one can end up with 
very specific, easy-to-evaluate expressions for the solution, useful, e.g., in the 
tsunami wave problem. 
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1. Introduction 

In the present paper, we use the theory of functions of noncommuting opera- 
tors pF-^ , aka noncommutative analysis (which can be viewed as a far-reaching gen- 
eralization of pseudodifferential operator calculus) , to solve an asymptotic problem 
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for a partial differential equation and show how, starting from general constructions 
and operator formulas that seem to be rather abstract from the viewpoint of dif- 
ferential equations, one can end up with very specific, easy-to-evaluate expressions 
for the solution, useful, e.g., in the tsunami wave problem. 

We consider the Cauchy problem with zero initial data for a 2D wave equation 
with variable velocity and with right-hand side localized near the origin in space and 
decaying in time. One physical interpretation of this problem is that it describes, 
in the linear approximation, the propagation of tsunami waves generated by local 
vertical displacements of the ocean bottom (see [4l-l9l ll7H19| and also |11H16| and 
the bibliography therein). Normally, the diameter of the region where these dis- 
placements occur (some tens to a hundred of kilometers) is much smaller than the 
distance traveled by the waves (thousands of kilometers), and their ratio, jj,, can 
serve as a small parameter. Accordingly, we are interested in the asymptotics of the 
solution as /I ^ 0. In the simplest piston model of tsunami generation, the bottom 
displacement occurs instantaneously at t = 0. This corresponds to a right-hand 
side of the form d'{t)v{x), where 6{t) is the Dirac delta function, and the problem 
is immediately equivalent, via Duhamel's principle, to the Cauchy problem for the 
homogeneous wave equation with initial data v(x) for the unknown function itself 
and zero initial data for its t-derivative. Fairly explicit asymptotic solution formulas 
suitable for easy implementation in Wolfram Mathematica |27j were constructed 
and analyzed for the latter problem in |11H16| on the basis of a generalization of 
Maslov's canonical operator jTl l20] . Now assume we wish to take into account the 
fact that the ocean bottom displacement evolves in time rather than happens in- 
stantaneously. Then it is natural to consider a right-hand side of the form g'{t)v{x), 
where g{t) is some smooth approximation to the delta function. An analysis shows 
that the solution can be represented as the sum of two parts, a propagating part, 
which travels along the characteristics, and a transient part, which is localized in 
the vicinity of the origin and decays in time. The propagating part can further 
be represented as the solution of the Cauchy problem for the homogeneous wave 
equation with initial data obtained from v{x) by application of certain functions 
/(L) of the spatial part L of the wave operator, where the corresponding symbols 
/(^) are given by simple formulas expressing them via the Fourier transform of g{t). 
These initial data, also localized near the origin, will be referred to as the equivalent 
source functions. The transient part of the solution is given by a formula similar to 
those for the equivalent source functions with the only difference that the function 
/(^) additionally depends on time as a parameter. The transient part is apparently 
not so important in tsunami wave analysis, but nevertheless it might be useful 
from the viewpoint of satellite registration of tsunami waves [17H19j . Since, as was 
mentioned above, the asymptotic formulas for the solution of the Cauchy problem 
with localized initial data for the homogeneous wave equation are already known 
from |llfll6) , we see that the only remaining thing is to compute f{L)v for all 
these functions /(^). It is here that noncommutative analysis comes fully into play. 
Note that L is an operator with variable coefficients, and so computing the function 
f{L) efficiently may prove quite a challenging task. However, all we actually need 
is the asymptotics of f{L)v, and methods of noncommutative analysis permit one 
to prove that f{L)v — /(Lo)w plus an asymptotically small remainder, where Lq is 
obtained from L by freezing the coefficients at the origin. Now computing f{Lo)v 
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is a breeze, because /(^o) is conjugate by the Fourier transform to the operator of 
multiphcation by the function /(crLo(p)), where CLgip) is the symbol of Lq. 

The one-dimensional counterpart of the problem studied in the present paper 
was considered in |21| . In the two-dimensional case, the results were announced 
in |22| . where the proofs were partly only sketched and partly absent altogether. 
Here we develop and refine these results and give complete proofs. Finally, note 
that we deal with the setting in which the wave propagation velocity is assumed to 
vanish nowhere. The case it which it vanishes (as it happens on the coastline in the 
tsunami run-up problem) is much more complicated. The asymptotics of solutions 
of such degenerate problems in some special cases was considered in the spirit of the 
approach of |llH16j in |23H25j (see also references therein); in the present paper, 
we restrict ourselves to wave propagation in open ocean. 

The outline of the paper is as follows. In Sec. [21 we give a detailed statement 
of the mathematical problem and write out well-known formulas expressing the 
solution in operator form. Using these formulas, we split the solution into the 
sum of the propagating and transient parts. Section |3] presents simple formulas for 
the asymptotics of the solution. The proofs of the theorems stated in this section 
depend on the results presented in Sec. SI which is the most important part of 
the paper and where the asymptotics of the equivalent source functions and the 
transient part of the solution are computed with the use of the noncommutative 
analysis machinery. Finally, Sec. [5] provides two simple examples; all computations 
and visualizations in these examples have been done with Wolfram Ma.thema.tica. 

2. Exact solution 

2.1. Statement of the problem. Consider the Cauchy problem for the wave 
equation 

with the initial conditions 

(2.2) ?7|t=o=0, 77t|t=o=0, 

where x = (xi, X2) G R^, r/ = ri{x, t) is the unknown function, c(x) is an everywhere 
positive smooth function stabilizing at infinityQ and the right-hand side Q = Q{x, t) 
depends on two parameters X, fi > and has the form 

(2.3) Q{x,t)^X'g',{Xt)v(^^^ 

with some smooth real functions V{y), y G R^, and goir), r £ [0, 00), such that 

(2.4) |X^(")(y)| <a(H-|y|)-l"l— , |a| = 0,1,2,... , 

(2.5) 5o(0)=0, / go{T)dT = l, |<?^'^(r)| < C^e-^^^ fc = 0,l,2,..., 

Jo 

for some >f > 1, > 0, and positive constants Cq and Cfe. 

Remark 2.1. One can also consider the case in which go{T) decays as some 
(sufficiently large) negative power of r as r — > cx). In this case, the estimates are 



That is, c{x) = const > for sufficiently large \x\. 
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somewhat more awkward, and we restrict ourselves to the case of the physicahy 
natural exponential decay (|2.5|) for the sake of clarity. 

Our aim is to find the asymptotics as /i of the solution of problem (|2.1|) 
on an arbitrary finite time interval uniformly with respect to A in the region 



This will be done in Sees. [3] and 21 and in the present section we deal with the exact 
solution of the problem. 

2.2. Physical interpretation and examples of right-hand sides. First, 
speaking in terms of the physical interpretation given in the introduction, let us 
explain the meaning of the parameters A and /i and condition (j2.6|) . The right-hand 
side Q{x, t) describes the time evolution (the factor X^gQ{Xt)) and the spatial shape 
(the factor V{x/iJ,)) of the perturbation (the tsunami source). In view of (|2.5|) . A 
characterizes the decay rate of the perturbation, so that 1/ A ~ to, where tg is the 
mean lifetime of the perturbation. The small parameter /i characterizes the source 
size tq, fj, ^ vq. We see that the product A/i = ro/io has the dimension of velocity 
and rewrite condition (j2.6|) in the form 

2.7 — = < WQ, 

Am ra 

where co = c(0), the wave propagation velocity at the origin, is taken to represent 
the typical wave propagation velocity in the problem and ujq is some dimensionless 
constant. This has a very clear meaning: the waves excited by the perturbation 
cannot travel too far before the perturbation dies out; they only cover a distance 
(co^o) of the same order of magnitude as the diameter rp of the perturbation region. 
We introduce the ratio 



so that condition ()2.7p (and hence p.6|) ) becomes 
(2.9) uj < ujQ. 

Mathematically, condition ()2.9p means that the parameter A is large (at least of 
the order of/i~^)asM— >-0. Note that, in view of the first two conditions in (|2.5p . 
XgoiXt) S{t) and X'^g'^iXt) -j> S'{t) as A -J> oo. 

In what follows, the dependence on the parameters A and /i is sometimes not 
immediately important to the argument, and in such cases we often "hide" these 
parameters by using the notation 



Next, let us give specific examples of right-hand sides Q{x,t). In practice, 
the actual ocean bottom displacement is known neither in much detail nor very 
precisely, because the corresponding measurements are impractical or impossible 
(cf. [17H19| V This results in certain freedom, which can be turned into an ad- 
vantage. Namely, when constructing the function Q{x,t) = g'{t)v{x) to be used 
in the analytical-numerical simulation according to the model (|2.1I) . one should 
take ansatzes that, on the one hand, fit the general information available about the 
source shape and evolution and, on the other hand, can be handled efhciently in the 
computations. (The latter includes the requirement that these functions, as well as 



(2.6) 



A/i > const > 0. 



(2.10) 
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Figure 1. The function V{y) with 61 = 1 and 62 = 4 rotated by 
the angle 6 = tt/IO (left) and its Fourier transform V{p) (right). 



their Fourier transforms, be given by closed-form expressions, which permits one to 
reduce the amount of numerical computations in favor of the less time-consuming 
analytical transformations.) 

A useful class of functions V{y) satisfying (|2.4|) is described by the expres- 
sion [T0l [T4l[T6] . generahzing [SjUli], 

\ 2 / \ 2x -3/2 

yi\ , / 2/2 ^ ^ 



^'■''^ ^(^)-n'ney ■ v^2 

where A, bi, and 62 are real parameters. The Fourier transform of this function is 
remarkably simple. 



(2.12) V{p) = A6i62e"V^!p!+lpi, 

One can further apply a differential operator 



P = P \ 

ydyi dy2 

with constant coefficients to the function V and then rotate the coordinate system 
by some angle 9, thus obtaining a broad variety of functions of the form 



Vp,e{v) = [PV]{T[e)y), T{0) 



cos 6 sin 9 
- sin 9 cos 9 

satisfying (|2.4I) . Such functions model elliptic-shaped sources of various eccentric- 
ity and various direction of axes with a wavy relief depending on the differential 
operator P (see ^ll, 15^). Figure [T] shows the graph of V{y) rotated by an angle of 
tt/IO and of its Fourier transform. 

Let us also give two examples of functions 50 (0 satisfying (|2.5p . 

(2.13) (a) 5o(r) =ae-"(sin(ar + 0o)-sin0o), (b) goir) = P{t), 

where a > and (j) are real parameters, a = (a^ -I- l)/{acos(f>o — a^sin^o) is a 
normalizing factor, and P{t) = X]fc=i(^0~^^fc''''^ is a polynomial of degree n with 
Efc=i'Pfc = l (see FigEl). 
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Figure 2. Examples of g{t) = Ago(At) for A = 1,2,3,4: go{T) = 
e~'^(sin(aT + ipd) — sin(po) (upper diagram); go{T) = e^'^(0.2r + 
0.4t^) (lower diagram). 



2.3. Operator solution formulas and energy estimates. We denote the 
spatial part of the wave operator in (j2.ip by L; thus. 



(2.14) 



Lu ~ — 



d 
dxi 



cHx) 



du 
dxi 



dx2 V 9x2 



-(V,c2(x)V)w. 



The operator (|2.14l) (with domain W|(R^)) is a nonnegative self-adjoint operator 
on L^(R^). Let D — ^/L be the positive square root of L. Problem p.ip becomes 

(2.15) r,"it) + D^r^it)^g'{t)v, ^ Vt\t=o ^ 0- 
Duhamel's formula represents the solution of (|2.15p as the integral 

(2.16) v{t)= f w{t,T)dT, 

Jo 

where w{t, r) is the solution of the problerrH 

(2.17) w[[{t,T)+D^w{t,T)^0, w\t^r=9{T)v, Wt\t=r=0. 

Indeed, the function (|2.16l) satisfies (I2.15|) . because 

/■* d 
r]"{t) + DMt)^ J {w['t{t,T)+DMt,r))dT + -{w{t,t)) +w't{t,t) 

^g'{t)v, 7,(0) =.0, 7y'(0)=g(0)z; = 
in view of (|2.5p . Now we can use the general solution formula (e.g., see |26[ p. 191]) 

(2.18) u{t) = cos{Dt)uo + D^^ 8in{Dt)ui 
for the abstract hyperbolic Cauchy problem 

(2.19) u"{t) + D^u{t)^0, uU^o^uo, utU^o^u, 
and write 



(2.20) 



coa{D{t — T))g{T) dr 



V = Re 



■^The standard Duhamel formula would give w\t=T = and wt\t=T = g'i^)^ in 112.171 1. but 
we have made use of the special form of the right-hand side. 
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Here we have used the fact that .g(r) is real- valued; the real part of an operator A 
is defined as usual by Re A = ^{A + A*). Formula (j2.20p is the desired abstract 
operator formula for the solution of problem (j2.ip . 

Remark 2.2. Since the operator D is self-adjoint, it follows that the expressions 
cos{Dt), s\Ti{Dt)/D, and e*^*, occurring in (|2.18p and (|2.20l) . are well defined in 
the framework of functional calculus for self- adjoint operators as functions f{D) 
with bounded continuous symbols /(^) = cos^t, f{£) = sin^i, and /(^) = e'^', 
respectively. Moreover, e'^* is none other than the strongly continuous group 
of unitary operators generated by D, Re(/(Z3)) = (Re/)(D), and, for "good" 
functions /(O; the operator f{D) can be defined not only via the integral over the 
spectral measure but also via the Fourier transform as 

/(i?)u=^(/(r),e-^w), ^eL2(R2), 
V zvr 

where / is the Fourier transform of / and the angle brackets stand for the value of 
the distribution fir) on the L^(R^)-valued function e"^u. 

Remark 2.3. The energy of the solution of the Cauchy problem (|2.19p is de- 
fined by the formula |261 p. 191] 

(2.21) E[u]{t) = \{\\u\t)f + \\Du{t)f) ^ \{\\u'{t)\f + iuit),Lu{t))) 

(where || • || stands for the norm and (•, •) for the inner product) and is 
conserved in the course of time. Hence, in view of ()2.10p and the estimates (|2.4p . 
the solution of problem ()2.17p (with t viewed as a parameter) satisfies 



X 



dxidx2 



^'^^^ ( c\^ly)\yV{y)\'dv^dy2^^gl{\T)cl\\yV\\\l + 0{^,)) 



/R2 2 

as /i — > 0, where cq = c(0). Now it follows from (|2.16p that, with some constant C, 
fM(0< {^V^N(i,r)dr| +]^\\w{t,t)f 

< C i^xj\goiXr)\ drj + ^gliXt) \\vf 

<c|y^ \9o{r)\drj + !^gl{Xt) \\V\\' 

= 0(1) + 0(Ai'A2e-2-^*) = 0(1) + 0(w-2g-2.At). 

i.e., the energy of the solution is uniformly bounded as /i — ^ for all t > e > 0. 
(However, it may have a "spike" of the order of a;~^ for t ^ 1/A; of course, this 
is only important if a; <C 1.) In other words, we have chosen a physically natural 
normalization of the right-hand side of our problem. 

Remark 2.4. For the inhomogeneous wave equation 

(2.22) u"{t) + D^u{t) = F{t), 
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one has the energy identity 

£[u]{t)^£[u]{0)+Re f {F{t),u' {t)) dr, 
Jo 

which implies the well-known estimate^ 

(2.23) \\u{t)l^, + \\u'{t)l < C{t){\\u{0)l^, + \\u'{0)l + sup ||F(r)|lJ, 

Te[o.t] 

where || • stands for the norm on the Sobolev space H'^ — W|(R^); in particular, 
II ■ llo ~ II ■ II • Of the estimates (|2.23l) . the most important for us is the one with 
s = (corresponding to the sum of the energy integral and the norm), in which 
the main estimates for the norms of remainders in asymptotic formulas will be 
obtained. However, occasionally our argument involves estimates with different s. 

2.4. Solution splitting into propagating and transient components. 

Let us further transform formula (|2.20p to reveal the structure of the solution and 
represent it in a form suitable for subsequent computations. We have 

/ e*^(*-^)5(T)dT = / e'^^'-^'>g{T)dT- / e'^^'-^'> g{T) dr 

(2.24) ^° ^" 

= e'"' e-'"^g{T)dT- e-'^^ g{T + t) dr. 
Jo Jo 

Let H{t) be the Heaviside step function {H{t) = 1 for r > and H{t) = for 

T < 0), and, for t > 0, let 

/•OO /"OO 

G{U)= e-'^^g{r + t)dr^ e-'^^ H{T)g{T + t) dr 

Jo J -OO 



be the Fourier transform of \/2TrH{T)g{T + 1) with respect to the variable r. (Note 
that G(^, 0) = -\/27r5(f ), where the function gij) is assumed to be extended by zero 
for the negative values of r.) Then formula (I2.24[) can be rewritten as 

/ e'-D(*-^)g(r)dT = e'^'G{D,{)) - G{D,t) = V2^e'^*g{D) ~ G{D,t), 
Jo 

and accordingly 

(2.25) ri{t) = riprop{t) + r]trans{t), 

where 

(2.26) r]prop{t) = V2^Re{e'°*g{D))v 

= V2ncos{Dt) Reg{D)v - V2nsm{Dt) lmg{D)v, 

(2.27) r]trans{t)^-Re{G{D,t))v. 

The function i]prop{i) given by (|2.26p is the solution of the Cauchy problem for 
the homogeneous wave equation 

(2.28) u"{t) + D'^u{t) = 
with the initial data 

(2.29) u^ = V2TTReg{D)v, ui =^ -V2Trl-aig{D)Dv. 



■^Their derivation takes into account the fact that the norm \\u\\ is equivalent to the norm 



II (1 + L)''/^n|| by virtue of the conditions imposed on the velocity c{x). 
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(This follows from the comparison of p.26p with (|2.18p .) Hence this function will 
be called the propagating component of the solution, and the initial data ()2.29|) for 
Vpropit) will be called the equivalent source functions. We shall see in Sec. I3.3l that. 
exactly as one should expect, rjprop{t) propagates along the characteristics. 

The function rjtrans (t) given by (j2.27p will be called the transient component of 
the solution, because it exponentially decays as At — ^ oo, as shown by the following 
proposition. (We shall also see in Sec. 13. ll that r]trans{t) always remains localized 
near the origin.) 

Proposition 2.5. As ji — > 0, the propagating component satisfies the estimates 

hp™p(0lll-O(l), Hrop(t)\\^0{l), 

and the transient component satisfies the estimates 

where v is the constant in condition (j2.5l) . 

Proof. We will estimate the transient part (|2.27p of the solution directly and 
the propagating part (I2.26P via the Cauchy data (|2.29p by using the energy es- 
timates. Formulas ()2.27p and ()2.29p involve the real and imaginary parts of the 
operator G{D,t) applied to the original right-hand side source function v. (Recall 
that g{D) is a special case of G{D,t) for t = 0.) Thus, we need to estimate the 
operator G{D,t). Note that, for an arbitrary bounded measurable function /(C), 
one has 

(2.30) ||/(i?):i/"^i/l <sup|/(0|, \\f{D):H'^H^\<Csnp\f{0\ 

with some constant C independent of /Q Thus, we need estimates for the function 
G{^,t). Since g{T) — XgoiXr), we have 

(2.31) G(C,t) = Go(e/A,At), 
where 

nOO 

(2.32) Go{^,t) ^ / e~'i^go{r + t)dT 

Jo 

is the Fourier transform of the function \/27r_ff (T)go(''' + t) with respect to t. By 
Lemma 12.61 below, we have 

1^2^5(01 -|Go(^/A,0)|< Coo 
and hence, by ((2:291) and (|230)) . 

Il^olli - V2^\\Re g{D)v\\, < CCqo M, = 0(1), 
||mi|| V2^\\lmg{D)Dv\\ < Coo \\Dv\\ < C \\v\\^ = 0(1), 

^Indeed, the first estimate is obvious, because the operator D is self-adjoint on = L2(R^)- 
To obtain the second estimate, we replace the norm on by the equivalent Hilbert norm = 
(n, {1 + L)u)^^'^ (cf. Remark l2.4l l: then the operator D becomes self-adjoint on H^, and the second 
estimate follows. 
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because v — V{x/fi) and hence = 0(1) (cf. the computation m Remark l2.3p 
Now the energy estimates (|2.23p for s = give the desired estimates for rjpropit) 
The estimates for the transient part go as fohows, again with the use of Lemma 

Wvtransim^ = 1 1 Rc G(-D, i)« 1 1 1 < C sup| Go (C/ A, At ) | 



<CCooe-''^'\\v\\,^Oie"'^'), 



\Vtransi^)\\ ~ 



dG, , 
Re—{D,t)v 



< sup 



dt \y 

< CoiXe-''^' \\v\\ = OifiXe-"^') = ©(w-^e"''^*), 
because v = V{x/fj,) and hence ||w|| = 0{^). This completes the proof. □ 

The foUowing lemma establishes the estimates for Go(^,t) used in the proof 
given above and also estimates that will be useful below. 



Lemma 2.6. The function Go{S,,t) satisfies the estimates 



(2.33) 



grn+kQ^ 



<Ck,ne~''*il + \m 



-k-1 



k,m = 0, 1,2, . 



dt"^d^'' 

with some constants Ckrn- For t — and m — 0, one has the better estimates 
d^Go, 



(2.34) 



-(e,o) 



<Cfco(l + ICI)""^ 



fc-2 



0,1,2,. 



Proof. First, let us prove the estimates ([05)1 and ([^31)1 for |^| < 1. Then 
we have 

^ I i-^rre-'^^^9t\r + t)dr 



< Cme 



by (|2.5p . whence the claim follows. Now let |^| > 1. Then we write 



T''e~''^dT 



O^Gc 
9t™ 



N poor iN 



dT^ 



gt\r + t)dr 



for some integer N > k + 1 and then integrate by parts TV times, thus obtaining 



(2.35) 



9" Go 



N 



-N 



1=1 



e~^^^g^^^''\T + t)dT. 



Next, we differentiate both sides of (12.351) k times with respect to ^, which gives 

Q^n+kQ^ N 



dt"^d^'' 



(i-i) 



k\ ( ? + ,s - 1)! 



r'=--e-^«^g^"+'^)(T + t)dT. 



Here all factors g^~^^ {t) and the integral are bounded in modulus by const • e""^* 
by virtue of (|2.5I) . and the smallest power of on the right-hand side is <^~'^"^, 
which implies the estimate (|2.33p . For t — Q and m = 0, the smallest power of ^"^ on 
the right-hand side is C"'^^^, since 50 (0) = 0, and we have the estimate (|2.34p . □ 
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3. Asymptotics of the solution 

In this section, we describe the asymptotics as /i — > of the solution 7/(t) = 
Vpropit) + rjtransit) of problem (|2.ip . (I2.2p . In all theorems in this section, we 
assume that all conditions stated in Sec. 12.11 are satisfied. Recall that the problem 
also contains the large parameter A, which is related to fi by the condition uj < ujq 
(see p.Qp ). where w — co(A/i)^^ (see ()2.8|) '). If u can be treated as a second small 
parameter (i.e., the distance traveled by the waves in the lifetime of the source is 
much smaller than the source diameter), then additional Taylor series expansions 
in uj lead to further simplifications in the asymptotic formulas. 

3.1. Asymptotics of the transient component. The asymptotics of the 
transient component rjtransit) of the solution as /i — > is given by the following 
theorem. 

Theorem 3.1. One has 

(3.1) Vtrans{x,t) ^ f f Re Go{uj\p\, Xt)V {p)e'P'' ^ dpi dp2 + R{t) , 



27r „ „ 

or, in the polar coordinates {r,ip), x ~ rn{ip), where n{ip) ~ (cosip,sinip), 

(3.2) rjtrans{rn{ip)) ^ / pReGo{ujp,Xt) 

27r Jo Jo 

X V{pn{i;))e'''P''°''^'''-^'>/^'dpdi; + R{t), 
where the remainder R{t) satisfies the estimates 

(3.3) \\R{t)\\^^0{p,e-''^'), \\R'{t)\\ =Oi^lio~^e~''^*), ^^0. 
Proof. Consider the operators 

(3.4) L(") - -clV\ = (i("))i/2. 

Thus, L^^'' is obtained by freezing the coefficients of the operator L at the origin, 
and D'^^^ is just the positive square root of the positive self-adjoint operator L'^^\ 

Lemma 3.2. One has 

(3.5) [ReG{D^"\t)-ReG{D,t)]v(^^^ ^ R{t), 

where R(t) satisfies the estimates p.3p . 

The proof of this lemma will be given in Sec. 01 Thus, the operator D in 
the expression (|2.27p for ritrans{x,t) can be replaced by the operator D^^^ with 
constant coefficients. Now it remains to compute ReG{D^^\t)V{x/ p). Since D'"' 
is an operator with constant coefficients and with symbol co|p|, it follows that 

(3.6) F(D(")) = J-ioF(cob|)o J- 

for any function -F(^), where J- is the Fourier transform and the middle factor on 
the right-hand side is the operator of multiplication by i^(co|p|). Thus we obtain 

G{D^''\t)v[j^=^ JjReGoiio\p\,Xt)V{p)e'P-/^dp^dp2 

(we have used the formula for G(^, t) and made obvious changes of variables), which 
proves Theorem 13. II □ 
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3.2. Asymptotics of the equivalent source functions. Now we proceed 
to the computation of the propagating component of the solution. It satisfies the 
Cauchy problem (I2.28|) . (|2.29p . and so a good starting point would be to compute 
the equivalent source functions (I2.29[) . Once we compute them (asymptotically) and 
prove that they are localized near the origin, we can use the methods developed in 
[11H16 ' to obtain the asymptotics of the propagating solution component. However, 
we would like to apply ready-to-use formulas from these papers rather than to write 
out new formulas based on the same ideas. The formulas in |llfll6) were obtained 
for the case in which ui, the initial data for the i-derivative of the solution, is zero. 
So we resort to the following trick. 

Proposition 3.3. The propagating solution component rjpropit) can be repre- 
sented in the form 

(3.7) Vprop{t) = Viit) + v'2{t), 

where r]i{t) and r]2{t) are the solutions of the Cauchy problems 

(3.8) 77i'(t) + = 0, rji\t=o ^ V2^RcgiD)v, Tyilt^o^O, 

(3.9) 7/^{t) + D^r]2 = 0, r]2\t=o^V2^D-^lmg{D)v, r/^|t=o = 0. 
Proof. The sum (|3.7|) obviously satisfies the wave equation (|2.28p . Next, 

?72|t=o = 0, (r]2)'|t=o = ??2li=o = ~D'^r]2\t=o = -V2TrDliiig{D)v, 

which shows that the initial conditions ()2.29|) are satisfied and hence completes the 
proof □ 

Thus, let us compute the asymptotics of the new equivalent source functions 

(3.10) 'nio^V2^Reg{D)v, i^a ^ D^^ lmg{D)v . 

Theorem 3.4. The equivalent source functions (|3.10|) have the following as- 
ymptotics as /i — )■ 0: 

(3.11) mo = Ui(^^^ +Ri, r?20 = C/2(^^^ +i?2, 

where the Fourier transforms of the functions Ui(y) and U2{y) are given by the 
formulas 

(3.12) UM = V2^Rego(^bl)^^(p), U2(j>) = V2^\-^ Im MMy(p) 

uj\p\ 

and the remainders satisfy the estimates 

(3.13) ||i?i||i = 0(M), \\R2\\2 = 0{^Ji). 

Proof. The proof goes along the same lines as that of Theorem 13. II Namely, 
we prove that the operator D in formulas p.lOp can asymptotically be replaced by 
D'^"^ and then compute the Fourier transforms of t/i and U2 using formula (|3.6I) . 
The latter computation is trivial, and we omit it altogether. As for the first part, 
it is given by the following lemma, which will be proved together with Lemma 13.21 
in Sec. IH 
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Lemma 3.5. One has 



v^[Reg(Z?)-Rcgp(°))]V^ - = Ri, 

(3.14) 



2tt[D-^ lmg{D) - (D^^^)-^ ImgpW)] (^^ j = R2, 
where Ri and R2 satisfy the estimates (|3.13p . 

This completes the proof of Theorem 13.41 □ 

Remark 3.6. If we replace rjio and 7/20 in the Cauchy problems for rji and ri2 by 
Ui{x/fi) and U2{x/fi), respectively, then the resulting error 5{t) in the computation 
of riprop{t) will satisfy the estimates 

Wsm.^oifi), \\6'it)\\=oi^i), ^i^o, 

uniformly on any finite time interval. Indeed, let us write S{t) — di{t) + 62(1), where 
Si and $2 are the errors in 771 and 772, respectively. Then, by virtue of the energy 
estimates (|2.23p . we have 

li^iwiii + \\sm\ < cit) \\6m\w = c{t) = o(m), 
ii^2Wiii + \m)\\ < c{t) \\5'm\ - c{t) \\dH2{q)\\ 

^C{t)\\D''R2\\<Ci{t)\\R2\\2 = 0{ii). 



Thus, the accuracy provided by Theorem 13.41 permits computing the propagating 
part of the solution modulo 0(/i) in the energy norm. 



3.3. Asymptotics of the propagating part. Remark 13.61 shows that, to 
compute the asymptotics of the propagating part of the solution of problem (|2.ip 
modulo 0(/i) in the energy norm, it suffices to solve problems p.Sp and p.9p 
asymptotically with the initial data replaced by the functions Ui{x/ fi) and U2{x/ fi) 
indicated in Theorem 13.41 Thus, we need to solve the problems 

(3.15) r,'l{t) + D^m=0, Tli\t=f)^Ui{x/n), 77il*=o=0, 

(3.16) 7^m + D^m=Q, m\t=o^U2{x/y), r,^|t=o=0. 

(We denote the new unknown functions by the same letters 771,2! this will not lead 
to a misunderstanding.) The initial data in these problems are localized near the 
origin, and hence the asymptotics of solutions of these problems modulo 0(/Lt) in 
all spaces H" can be obtained with the use of the approach developed in [llf 
114] and based on the Maslov canonical operator [ll [20j . Let us briefly recall this 
construction. 

3.3.1. Bicharacteristics, canonical operator and solution formulas. In the phase 
space p with the coordinates (x,p) = (xi, a:2,pi,P2), consider the Hamiltonian 
system 

corresponding to the Hamiltonian function % = |p|c(a;). This system determines the 
Hamiltonian phase flow g\^. Let n(7/;) = * (cos 7/', sin t/;). Consider the Lagrangian 
manifold = {p — n(i/)), x = aw{tp)}, isomorphic to the two-dimensional cylinder, 
where t/j S [0, 27r) and a G R are coordinates on Aq. By shifting this manifold along 
the flow (7^, we obtain the family of Lagrangian manifolds At — q^Aq, each of which 
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is equipped with the same coordinate system (tp^ a) as Aq. We take the point with 
coordinates {ip, a) — (0, 0) for the distinguished point on Aq and construct the 
Maslov canonical operator K'j^ W)20\ on each of the manifolds Af. (Here h 
is the small parameter occurring in the construction of the canonical operator; all 
Jacobians in the definition of are taken with respect to the coordinates a).) 

It follows from the results in |11H14I that the asymptotics of the solutions 
111,2 of problems (|3.15|) and p.l6p can be obtained as follows. Using the Fourier 
transforms (|3.12p of the equivalent source functions computed in Theorem 13.41 we 
introduce the following two smooth functions on Af, independent of t and a but 
depending on the coordinate ip and an additional parameter p: 



fiii^^p) ^ Ui{pn{tp)) = V2TTRego{ujp)V{pn{'ip)), 



Mi^,p) - C/2(pn(^)) = V2^A-ilm^^F(pn(V')). 

ojp 

Then the formulas in |111I13| give 

(3.18) mAt) = y^Re(^e-"/*^°°ifl/^(VP¥'i.2(V',p))dp) +0{p). 

Let us find the derivative rj2{t). By the commutation formula [20] for the canonical 
operator, we have 



^Re(e-"/^ K'i^/;[^'j-'H\^^^p^2{^,p)\ dp)+0{f,). 



27. 

But the Hamiltonian H is preserved along the trajectories of the Hamiltonian sys- 
tem, and hence "Hj^ — "Hj^ = c{an{ip)). It was shown in [11 that, modulo 
lower-order terms, one can set a = in the functions on Af. Taking into account 
the definition of ip2, we obtain 

77^(0 = Re (^e-T^ K^^/[^{-tlmgo{u:p)V{pn{i^)))] d + 0{p). 

Finally, we use the formula rjpropit) = rjiit) + ?72(i) and arrive at the following 
theorem. 

Theorem 3.7. The propagating part of the solution has the following asymp- 
totics: 

(3.19) Vpropit) = ^Re(^e~'f K^//[y^To {^p)V{pn(i;))] d p^ + R{t), 

where the bar stands for complex conjugation and the remainder satisfies the esti- 
mates 

(3.20) \\Rm,^Oip), \\R'm=0{^,) 



uniformly on any finite interval of time t. 
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3.3.2. Asymptotics near the front. Now let us compute the propagating part 
()3.19|) of the solution in more explicit terms. To this end, we need some geometry. 
Let {P{t,ip), X{t,ip)), ip e [0.27r), be the family of solutions of the Hamiltonian 
system (|3.17p with the initial conditions 

(3.21) p\t^o = n(V^), x\t^o = 0. 

For each t, the equations p — P{t,ip), x — X{t,ip), ip E [0.2tt), define a smooth 
closed curve Tt in the four-dimensional phase space Ri^^; this curve is called the 
wave front in R!^ p. The projection jt = {x = X{t, ip): E [0.27r)} of Vt into is 
called the front in the configuration space. In contrast to Tt , the curve 7t may well 
be nonsmooth; namely, it may have turning (or focal) points (in this case, = 
for some ■0) and points of self-intersection. Moreover, the front 70 at the initial 
time i = is just the point x ~Q. 

For each t, the function (|3.19p is localized in a neighborhood of the front 74 
[11H16| . Formula p.l9p provides the global asymptotics of the propagating part of 
the solution; i.e., this formula holds both near regular and near focal points of the 
front. The formula can be simplified in a neighborhood of any point of the front, 
but the simplified expression depends on whether the point is regular or focal. Here 
we restrict ourselves to the case of a neighborhood of a regular point. 

Take some time t and angle and assume that the point X{t,ip'^) G 7* is 
not focal; i.e., X^(i,'0°) ^ 0. In some neighborhood of X(i,-0°), we can introduce 
the local coordinates {ip,y), where y = y{x,t) is the (signed) distance between the 
point X and the front and V' = 4'{xj t) is determined by the condition that the vector 
X — X{t, ■ip{x, t)) is orthogonal to the vector tangent to the wave front at the point 
X(t,ilj{x,t)); in other words 

{x - X{t,ij{x,t)),X^{t,iP{x,t))) ^0. 

Set 

S{x,t) = {P{t,yj{x,t)),x-X{t,ij{x,t))). 

Next, we introduce the Morse index m{t,ip'^) of the trajectory X(t, i/''^), t <E (0,i], 
which is the number of zeros of the function \X^{t,iI)^)\ on the half-open interval 

re (0,t] m- 

It may happen that some region of points x where we intend to write out 
the asymptotics simultaneously belongs to several neighborhoods of the above- 
mentioned type, where the corresponding points X{t, lie on several (but finitely 
many!) distinct arcs of the front 74. (For example, this is the case if we study the 
asymptotics near a point of self- intersection of the front 74.) Then all these arcs 
contribute to the asymptotics at such points x, and we use an additional subscript 
j to distinguish these neighborhoods as well as all associated objects (■0°, ^{x,t), 
S{x,t), Morse index, etc.). Now from the results in |lllll2l[T5lll6| we obtain the 
following theorem. 

Theorem 3.8. In a neighborhood of the front jt but outside a neighborhood 
of the focal points, the asymptotic formula p.l9p for the propagating part of the 
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solution can be rewritten in the form 



(3.22) r]prop(t) = v^Re^ 



-g-i7rm(V°,t)/2 r 




j 



_^\X^{^,t)\ \lc{X{ij,t)) 



+ R{t), 



where 



■oo 



(3.23) 



F(z,^) 



— e 



iir/4 



R{t) satisfies the estimate (|3.20p . and the sum with respect to j is taken over all 
distinct arcs of jt contributing to the asymptotics at x0 

Remark 3.9. The factor 



includes the two-dimensional analog of the so-called Green law and the trajectory 
divergence related to the velocity c{x) (with height c'^{x) describing the bottom 
topography). The function F depends on the time and space shape of the source 
generating the waves |11H16] . Formulas (|3.19p . (|3.22p . and (|3.23p apply to any 
localized perturbation. 

4. Obtaining asymptotic expansions by noncommutative analysis 

The aim of this section is to prove Lemmas 13. 2l and r3. 51 Vaguely speaking, these 
lemmas state that the replacement of the operator D by the operator D'^"^ p.4|) 
with constant coefficients in certain expressions results in an 0{fi) error. However, 
it is much easier to deal with functions of the differential operators L and L*^"-' than 
with functions of their square roots, the pseudodifferential operators D — \/L and 
D^°^ = VU^. Hence in Sec. 14.11 we represent the latter functions via the former 
and accordingly restate the lemmas. In Sec. 14. 2[ we make all noncommutative 
computations. 

4.1. Eliminating the square roots. 

Lemma 4.1. The functions ReGo(C,i) and ^^^ImGo(^,i) are smooth even 
functions of ^ and hence smooth functions of . 

Proof. The function goir) is real-valued, and Go(C,i) = Go(-^,t) by (|2.32p . 
Thus, ReGo(-C,i) = RcGo(C,t) and ImGo(-^,i) = -ImGo(C,0; i-e-, the real 
part of Go is an even function of ^, and the imaginary part of Go is an odd function 
of ^. Hence the desired claim follows. □ 

Now let us introduce the functions 



More formally, for example, fix an e > 0; the intersection of 74 with the e-neighborhood of 
X can be covered by finitely many arcs of length < e; take the contributions of all these arcs. 




(4.1) 



/i(0 -RcGo(C^/^0), /2(0 =r'/'lmGo(e^/2^0) 
/3(e,i) =ReGo(?l/^^). 
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By Lemma 14.11 these functions are smooth for all ^, including ^ = 0. Formulas 
p.lO|) for the equivalent sources and (|2.27p for the transient solution component 
can now be rewritten as 



(4.2) 

Vtrans{t)^~f3{X-^L,Xt)v(j]. 

Indeed, for example, 

X-^h{X-^L) = A-i(A-2L)-i/2iniGo((A-2i)i/2,o) 

= D-HmGo{X-^D,0) = Im G'(L», 0) = V2^D-^ lmg{D). 



The following theorem is an equivalent restatement of Lemmas 13.21 and l3.5l in terms 
of functions of L and L^^\ (We write i?3(t) = ^Rit) to unify the notation.) 

Theorem 4.2. One has 

A(A-2l)f(^^) =/i(A-2L(0))y(^^) 



(4.3) x-'h{x-H)v^:!-j = x-\h{x-^L("^)v^j^j 

h{X-H,Xt)v(^^^ =/3(A-2L(o),Ai)^(^) +R3it), 
where the remainders satisfy the estimate 

PiIIi-O(m), ||i?2|l2=0(M), 

imrn, - 0(A.e-^^*), = 0(A.^Ae-^^*). 



The proof will be given below in Sec. 
We need some estimates for the symbols (14.11) . These are provided by the 
following lemma. 



Lemma 4.3. The following estimates hold for the functions (|4.ip : 



(4.5) 



dx^'dt" 



<Cfc„e--*(l + lcl)-l/'^^ 



fc = 0, 1, 2, . . . , where the Ckm cire some constants {in general, different from those 
introduced earlier). 

Proof. For k = 0, the desired estimates (|4.5p readily follow from (|2.33p and 
()2.34p : it suffices to replace ^ by ^^/^ (and use the fact that the functions fj given 
by (|4.ip are smooth and in particular continuous at f = 0). Next, note that if 
/(^) = F(^^/^), where F(C) is a smooth even function, then /'(^) — 5'(^^/^), where 
^'(C) = ^-F"(C)/C is again a smooth even function. Thus, it suffices to prove that if 
a smooth even function F satisfies estimates of the form 

|i^('^(C)l <'^fc(l + ICI)"'"'°, A = 0,l,2,... , 
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for some fco, then ^ satisfies the same estimates but with ko increased by 2 and 
with new constants dk, each of which is a finite hnear combination of the old ones. 
This is trivial for |C| > 1, and in the region |C| < 1 one can use the identity 

C'F' io = c'{F' iC) - F'{0)) = F"{9C) d9. a 

Jo 

4.2. Computation of the transient part and the equivalent sources. 

Now we will prove Lemmas 13. 21 and l 3.5l bv proving the equivalent Theorem 14. 2 1 Let 
/(^) be any of the functions /i(^), /2(0j ^^^d f^i^^t) given by (|4.ip or the function 
f4{xi,t) = df3{£_,t)/dt. We need to compute the difference 

(4.6) n ^ n{x, A, m) = (/(A-^L) - /(a-2l(")))f (^^) 

and estimate it in an appropriate norm. Let us make the change of variables x = fiy. 
In the new variables, (14.61) becomes 



(4.7) 7^=(/(L,)-/(L(")))y(y), 
where 

and uj = co/(A/i) is bounded by condition p.6p . 

To compute this difference, we use the machinery of noncommutative analysis. 
We refer the reader to [UlS] for details concerning the definition and properties of 

1 n 

functions of noncommuting operators and only recall that a function F{Ai, . . . , An) 
of (possibly, noncommuting) operators Ai, . . . An can be defined as follows in the 
particular case where the Aj are the generators of uniformly bounded strongly 
continuous one-parameter operator groups e'^^*, i S R, on a Hilbert space H: 

In 1 f ~ A A 

F{A^,...,An)u= / F{hM---.tn)e'^-'" ■■■e'^^'^udh---dtn, 

u G H, where F is the Fourier transform of the symbol F, which is assumed to 
satisfy certain conditions (e.g., see [2]) guaranteeing that the integral on the right- 
hand side is well defined. The numbers (Feynman indices) over operators indicate 
the order of their action: of any two operators, the operator with the smaller 
Feynman index stands to the right of the operators with the larger Feynman index 
in products. 

It follows by the zero-order Newton formula of noncommutative analysis (see [1] 
and m Theorem 1.8]) that 

(4.8) fiLy) f{Lf) = |(4,4°^)^Z~4^. |(4,4"))r, 
where 

-jT(l;i,t.2) T T 

04 ?1 - t.2 

is the first difference quotient of / and 
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Here we have denoted 

^0 

this function is uniformly bounded together with ah of its derivatives for z e R^, 
and 0(0) = 0. 

We further transform the right-hand side of (|4.8p as follows. 
Proposition 4.4. 

(4.9) ^(4,4°))T=|(4,4^)T+g(4,4°\4''^)[r,4°^], 
where 



772 (Cl''?2,6) t ^ 



3 



is the second difference quotient of f and 
is the commutator of T and L^'^ . 



Proof. The proof of ()4.9p mimics the derivation of the general commutation 
formula 2j Proposition 1.3] 

(4.10) 1/(B)-A/(I) = ^(B,B)[I^ 

of noncommutative analysis, with T playing the role of A and Ly"-* playing the role 
of B. Recall this derivation (e.g., see |2] pp. 52-53]). We need to compute 

[A, f{B)] ^ Af{B) - f{B)A - Af{B) - Af{B). 

The Feynman indices can be chosen independently for either summand on the right, 
and we can write 

[A, f{B)] - Af{B) - Af{B) = A{f{B) - f{B)) = MB - B)j,{B,B)- 

(Here we have used the identity f{x) — J{£) = {x — ^)||-(a;, ^), which is in fact the 
definition of Sf /S^.) Next, we move apart the Feynman indices over the S's, thus 
obtaining 

21 3Afl3 21 3Af04 ^ Xf 4 ^ Xf 1 3 

A{B - B)^{B, B)^A{B - B)^{B, B) = [A, B]^{B, B) = [A, B]f^{B, B). 



(In the middle, we have written A{B — B) — [A, B] using the fact that no other 
operators in the expression have Feynman indices in the interval [1,3].) Thus, we 
arrive at the desired commutation formula (|4.10p . 

1 

The derivation of (|4.9p differs from this only in that now, instead of f{B), 

A f 3 (f)\ 

we have -^{Ly,Ly ); i.e., there is an additional operator argument, Ly, but this 
argument does not invalidate the computation, because its Feynman number does 
not lie between those of A = T and B — L^'^ . □ 

Let us evaluate the commutator [T, 
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Proposition 4.5. One has 

[T, 4°)] - ^Ti, where \\T^ : i?^(R2) i?''-'(R^)|| < 
for all s with some constants independent of fj, as ^ — > 0. 
Proof. We have 



2 

.4 ' 



and it remains to recall that 0(z) is uniformly bounded together with all derivatives. 

□ 



By Propositions 14.41 and 14.51 we can write 
Accordingly, 

(4.11) n = {f{Ly) - /(40)))F ^AW + ^iBV, 
where 

(4.12) W = TV, A^^-l{Ly,Lf)., B =^^{LyJ^^\L^^^)T,. 

Let us estimate the expression (|4.1ip for f = fj, j = 1,2,3, 4. 
Proposition 4.6. One has V e H^CRy) for every s. 

Proof. This follows from the estimates (12.41). □ 



Proposition 4.7. For every s, one has W £ i7'*(R^) and 
Proof. We have 

(j){fiy) = fi{F{^y),y), 

where the vector function 

F{z)^ f ^{ez)de 

Jo 9z 

is bounded together with all derivatives, and hence for the function W = TV we 
obtain 

H,,,.^^(A(.,,„,,.,i5M,|.(,,,„,,,,2M). 

Since, by virtue of the estimates (|2.4p . the function yjdV{y)/ dyu lies in ^^^(R^) 
for every s, we arrive at the desired assertion. □ 

Proposition 4.8. Let f = fj, j = 1,2,3,4. Then for each s e R there exists 
a constant Cs independent of fj, ^ such that 

\\A: H^Rl) ^ H^RDW < Cs, \\B: H%R.l) ^ ff^-3(R^)|| < 
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for J ^1,2, 



\A: H'iTiiJ H'iKfJW < CsC-"', \\B: H^Rf^) H'-^(Rl)\\ < C -^"^ 



e 



for j = 3,4. 

Proof. We make use of the following representation of the /cth difference quo- 
tient: 

/ f^''\e^i, + --- + eu+i£.k+i)d6^---dek 



where J^*^) (p) is the Fourier transform of the fcth derivative (^) and 

Afe = {(01, . . . , Ok+i e R'-^+i : 01 + • • • + Ok+i - 1, 0, > 0, j = 1, . . . , fc + 1} 
is the standard fc-simplex. Hence 



dp d0i 



~ V2^ J 

Let us estimate the operators (I4.13p . To this end, we use the following lemma. 



/oo \ 
/" (p) e'P"' e'po^^y^ Ti e'^^^^y'' dp ) d^i dflj . 
-oo / 



Lemma 4.9. For each s, there exists a constant Cg independent of fi ^ such 

that 

||e'*^«: i7"(R2)->7J«(R2)|| <C„ e'^^y^ : H'{Rl) ^ H^RD < Cs 

for all t eR. 

Proof. For s = 0, the claim is obvious, because the operators L^,"-* and Ly 
are self-adjoint in L^(Rj,). For other values of s, one equips H''{Ry) with the 
equivalent norm || (1 -|- Lj,)"/^!*]], so that the operator Ly becomes self-adjoint. This 
norm depends on the parameter /z, but it is not hard to prove (for positive integer 
s by a straightforward computation, and for other s by duality and interpolation) 
that the constants in the inequalities specifying the equivalence of norms remain 
bounded as — >■ 0. The argument for L^"' is simpler, because the parameter /i is 
not involved. The proof of Lemma 14.91 is complete. □ 

Now we can finish the proof of Proposition 14.81 If / = /i,/2,/3, or f^, then 
it follows from Lemma 14.31 that the Fourier transforms of /' and /" belong to 
L^{R), and in the case of /a and f^ the i^-norm decays as e""^*. By combining 
this with Lemma [4.91 and with the estimate for Ti in Proposition 14.51 we arrive at 
the assertion of Proposition l478l □ 

By applying Propositions 14.61 14.71 and to formulas (14. lip and (|4.12p , we 
find that n = 0(/i) in aU iJ^(R2) for / = /i and / = /a and 7^ ^ 0(^6-"*) in ah 
H^'iRy) for f — fs and f — fi- Let us finally estimate the remainders Rj in (|4.3p . 
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We should take into account the additional factor A ^ for j = 2 and pass from the 
variables y to the original variables x — fiy. Since 

(4.14) ll"IL = ll"lli/.(R2) < M^"' ll"ll//.(R2) for ^ < 1 and s > 0, 

we arrive at the desired estimates (|4.4p . For example, for R2 we obtain 

Co 

(where the factor A^^ comes from (j4.2p and the factor /i^^^ — fi^^ from (I4.14p for 
s = 2). The estimates for Ri and R3 are similar. The proof of Theorem 14.21 and 
hence of Lemmas 13.21 and 13.51 is complete. □ 



5. Examples 

In conclusion, let us present two simple examples in which the asymptotics of 
the solution of the Cauchy problem (|2.ip , (|2.2p with a special right-hand side will be 
demonstrated. Namely, we use the right-hand side Q{x, t) = g'f){Xt)V (x / n) , 

where V{y) = A{1 + {yi/biY + (y2/&2)^)~^''^ is the simplest spatial shape fac- 
tor (j2.1ip and the function (7o(t) is given by one of formulas (a) (a sine source) and 
(b) (a polynomial source) in Eq. (j2.13p . 

Recall that the asymptotics of the solution is given by Theorem 13. 11 Eqs. (|3.ip 
and (|3.2p (the transient solution component) and by Theorem 13.81 Eq. p.22p (the 
propagating solution component away from the focal points). The transient com- 
ponent ritrans{x,t) and the wave profile F{z,ip) (see p.23p ) of the propagating 
component depend only on the right-hand side and on the parameters A,/i, and 
w; they are represented by integrals which, for our choice of the right-hand side, 
can be evaluated (or, in the case of the transient component, considerably simpli- 
fied) analytically. The other ingredients of the asymptotic formula p.22p for the 
propagating component (the phase functions Sj{x,t), the Lagrangian coordinates 
ipj{x, t), the Morse index m{ip'j, t), and the factors responsible for the Green law and 
for the trajectory divergence) depend on the solution of the Cauchy problem p.2ip 
for the Hamiltonian system p.l7p . which, except for the simplest cases, should be 
solved numerically. 

Accordingly, our exposition in both examples is as follows. First, we find the 
function Go{£,,t) (|2.32p . which plays a crucial role in all subsequent calculations. 
Then we write out the wave profile F{z,ip) and finally present the expression for 
the transient component rjtrans{x, t) of the solution. In the second example, we also 
numerically compute the trajectories and display snapshots of the solution obtained 
with the use of Wolfram Mathematica. 

The calculations are mostly carried out in polar coordinates, so let us rewrite 
formula (j2.12l) for the Fourier transform of V in the polar coordinates (p, ip), where 
p = /fn(-0) with n('(/;) = (cosV', sin?/j): 

(5.1) V{pn{'iP)) = Abib2e~P^'^'''\ where l3{^) = bj cos^ + sm^ip. 

5.1. The case of a sine source. Let 

goir) = ae"''(sin(aT + (po) - sin^o) 
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where a — (a^ + l)/{a cos (j>o — o? sin 0o) is a normalizing factor. By evaluating the 
integral in (|2.32|) . we obtain 



We see that Go(^, i) is a rational function of ^. Moreover, a routine computation 
(which we omit) shows that it can be represented in the form 

(5.3) Go(e,t) =^<Z™(i)(^m(C')+«eQm(e')), 



where -Rm(C) ^-nd Qm{C) ^-re rational functions with real coefficients and with de- 
nominators nonvanishing for C ^ 0. This is, of course, consistent with the assertion 
in Lemma 14.11 concerning the parity of the real and imaginary parts of Gq. As to 
50 (C): have 

To evaluate the wave profile F{z,ijj) of the propagating solution component, 
we substitute the functions (|5.ip and (|5.4p into formula p.23p and obtain 

^, aAbib2e-"/^ 



X 



27rcj3/2 

V 1 + ia — i|0 1 — m — ifl 1 — ipJ 







27rcj3/2 



+ ia — ip 1 — ia — ip 1 — ipj 
'-e-'^no{uj-\l3{4>)-iz),l + ia) 



- -e''^°Io(c^-H/3W - ^z), l-ia)- lo (c^-i(/3(^) - iz), l) sincp^ 
where the integral 

(5.5) Io(Ci,C2) = / ^ Ci,C2 e C, RcCi > 0, argC2 ^ -, 

Jo C2-ip 2 

can be expressed via the complementary error function 

erfc(w) = — ^ / e"" cZi; 

by the formula 

Io(Ci, C2) = ^ + e-'^^/Vy^e'^i^^ erfc (e^/^^^). 
vCi 
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To evaluate the transient term of the solution, we substitute the functions ()5.ip 
and (15.21) into (13.21) and obtain 



r]trans{rn{ip)) 



27rw2 

a — i 

a + i 



Re 

,-i(at+0o) 
,i{at+<po) 



aAbib2e 
27rcj2 



-At /•27r 



Re 



— (sin(aAt + (/iq) ~ sin(/)o) + sin^o 
z 

p + a~ i 
p- a- i 



P-i 



dip + 0{pl) 



(sin(aAt + <t>o) ^ snKpo) 



+ sin 00 26"'^ (tt + 2iCi(z) - 2Si(z)) 

+ ^_l^-^{'-t+M^{—'^)^Ei{{a - i)z) 

where z = z(r, 93, ■(/;) = '^^"'^ (/^(V') ~ irfi~'^ cos(V' — f/?)) , Re(z) > 0, and 



n+oo ~t 

Ei{z) EE I —dt, Ci(z) 



cos t 



dt, Si (z) = / dt 



sint 



5.2. The case of a polynomial source. Now let 

where 

n 



fc=i 



is a polynomial of degree n with coefficients Pk such that Po = and X]fc=i -Pfc = 1- 
Let us use formula (|2.32l) for Go(C,t). Since 



{t + r)*^ dr = e-'\t + i 



1 



5 



it follows that 



(5.6) Go{i,t) = e-'P[t + i^ 



1 



30(0 = 



1 



27r V^C/l + ie' 



1 



and we see that Go(^,i) again has the form (|5.3p . Using ()5.ip . (|5.6p and p.23p . we 
evaluate the wave profile of the propagating part of the solution as follows: 



2^tj3/2 



C2 = l 



.Ablb2\/TT 



,iCi 



where Io(Ci,C2) is the integral (|5.5 
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Remark 5.1. In both examples, one can prove that the following asymptotic 
formulas hold for the functions F{z,'>p) for small uj: 

This means that for small uj the solution of the inhomogeneous problem (corre- 
sponding to "sources stretched in time" ) passes into the solution of the homoge- 
neous problem (corresponding to "instantaneous sources"). 

Let us compute the transient term of the solution for the case in which P(t) is 
a second-order polynomial; then 



t f P2ty2 + {Pi - P2)t - P, 2P2t + 2Pi-3P2 , 4P2 



— i^e 

For the transient term, we find 

_ _\2 
'Utrans AC 



i+e {i+e)^ i^+e)^ 

JP2t'/2 + Pit 2P2t + 2Pi-P2 , 4P2 



2„-At ip^xhy2 + {Pi-P2)\t-pAei"' 



{2P2\H+{2Pi - 3P2)A2)e2(-) +4P2A4e3(-) 



where 



27rA^'^ V (1 + [^^mrr 

If we pass to the polar coordinates by setting y — m{Lp) and p — pniip), then we 
obtain 

Abib2 f°° f^^ pg-pifiW—ir cos(i^-^)) 

e.(m(c,),^) = ^y^ ^T^^V)I dpdi,. 

Here one can evaluate the integral over p. For k = 1,2,3, we obtain 

ei{rn{ip),p) ^ ^^^2^2 J ci^(-cos(z)Ci(z) + isin(z)(7r-2Si(z))), 

e2(rn((^),/i) = dv(2-2zsin(z)Ci(z)-zcos(z)(7r-2Si(z))), 

e3(rn(v9),/i) = / dv/4-zsin(z)(^z + 2Ci(z)-2zSi(z)) 

+ z cos(z)(-7r + 2zCi(z) + 2Si(z))^ , 

where z{il}) = u:~^{P ~ ircos{ip — (f)). 

An illustration of the solution given by the sum of propagating and transient 
terms in the second example is shown in Fig. |31 Here the propagating part is 
calculated for the constant velocity c{x) = cq = 1, and other constants are hi = 
1,62 = 2, A = 1,/i = 0.1, Pi = 0, P2 = 1. The first four snapshots are taken 
at small times t = 0.3,0.7,1.0,1.5 to show how the transient term behaves, and 
the last three snapshots are taken at large times t — 1.5,4.0,6.5. At i = 6.5, 
the transient term practically disappears, while the propagating part continues its 
motion. The function 170 and the wave profile for Pi — —2, P2 — 3, and various A are 
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■t= 0.3 ; Pi= ; P2= 1 ; O.U /t^ 0.1 ; Pi= ; Pj^ 1 ; O.U 

1 ; 2 ; A= 1 ; A= 5 J Ui= 1 ; 2 ; A= 1 ; A= 5 J 




■t= S.5 , Pi= ; P2= 1 , fi= 0.1 \ 
.Sl= 1 , ^2= a ; A= 1 , A= 5 J 




Figure 3. Sum of waves jyprop + ?7t: 



compared in Fig. 21 For small A, the wave profile has the form that "reproduces" 
the shape of the function go, while for large A the wave profile is almost the same 
as for go — 5{t). 
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g(t) 

20 r 

15 ^1 1 

10 r| I 



5 ? 1 




-5li 

(T= 1 , 3 ; Pi= -2 ; P2= 3 ; //= 0.1 ; A = Afi = 0.5, 1, 5, 100) 



'/prop 
1.2 r 



1.0 - II „ 

w 
I I.' 



I ) 




Figure 4. Examples of profiles of propagating waves. 
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